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M. 
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V 


L 
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reference  configuration 
vectors  in  1R3  based  at  X  €  8 
deformation 

infinitesimal  deformation 
strain 

all  deformations 
deformation  gradient 
transpose  of  F 
Cauchy  Green  tensor 
Stored  energy  function 
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all  loads  with  total  force  zero 
all  linear  maps  of  Tx3  to  1R'5 
symmetric  linear  maps  of  TXB  to  Tx5  . 
(Q  €  lOR'5^'5)  !  QTQ  =  I,  det  Q  =  1} 
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symmetric  elements  of  M, 

a 

skew  symmetric  elements  of  M„ 

.■) 

infinitesimal  rotation  about  the  axis  v 


equilibrated  loads 


astatic  load  mat) 


k  :  L  -*■  M. 
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sym 
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astatic  load  for  a  load  2, 

non  singular  part  of  k 

skew  viewed  in  load  space 

sym  viewed  in  load  space 

$(<f0  =  (-DIV  P,  P-N) 

{u:  3  -*-]R3|u(0)  =  0,  7u(0)  £  sym} 

image  of  C  near  I„  under  t> 

B  sym  B 

N  is  the  graph  of  F 
F(\,Z)  =  F(X2)/A: 
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§1.  \lntroduction 
^ - 

The  purpose  of  this  paper  is  to  study  the  traction  problem  in 
three  dimensional  nonlinear  elasticity  using  geometric  techniques  and 
singularity  theory.  The  first  two  oapers  in  the  series  deal  with  dead 
loads  and  configurations  that  are  nearly  stress  free.  As  was  shown  by 
Signorini  [1930]  and  Stoppelli  [1958],  this  problem  has  nontrivial 
solutions.  However,  their  analysis  is  incomplete  for  three  reasons. 
First,  their  load  was  varied  only  by  a  scalar  factor;  in  a  full 
neighborhood  in  load  space  of  a  load  which  has  an  axis  of  equilibrium 
there  are  additional  solutions  missed  by  their  analysis.  Second, 
their  analysis  is  only  local  in  the  rotation  group,  so  additional 
nearly  stress  free  solutions  are  missed  by  restricting  to  rotations 
near  the  identity.  Third,  some  classes  of  loads  with  a  degenerate  axis 
of  equilibrium  were  not  considered.  This  paper  completes  their  analysis 
by  treating  these  questions  as  well  as  stability.  The  complexity  of 
the  answer  is  indicated  by  the  fact  that  near  certain  types  of  loads, 
we  find  up  to  40  distinct  solutions  that  are  nearly  stress  free.  Our 
constitutive  hypotheses  on  the  stress  tensor  are  'generic';  for  a 
degenerate  stress  tensor  there  can  be  even  more  solutions. 

The  literature  on  this  problem  is  very  extensive,  going  back 
to  Signorini  in  the  1930's.  Our  primary  sources  have  been  Stoppelli 
[1958],  Grioli  [1962],  Truesdell  and  Noll  [1965],  Van  Buren  [1968], 

Wang  and  Truesdell  [19"],  and  Capric  and  Podio  Guidugli  .1974], 

However,  none  of  the  references  beyond  Stoppelli  [1958]  proves  any 
of  the  theorems  dealing  with  nontrivial  cases;  i.e.  loads  with  axes 
of  equilibrium.  However,  Grioli  'I960]  is  a  convenient  reference 


for  the  statements. 


The  outline  of  this  first  part  of  the  paper  is  as  follows.  Our 
notation  for  nonlinear  elasticity  and  the  problem  near  a  natural  state 
is  formulated  in  Section  2.  In  Section  3  the  basic  properties  of  the 
astatic  load  are  reviewed  and  developed.  The  problem  is  reformulated 
with  special  reference  to  the  global  aspects  of  the  rotation  group  in 
Section  4  and  introduces  the  bifurcation  equation  which  plays  a  crucial 
role  throughout  the  paper.  Section  S  treats  loads  with  no  axis  of 
equilibrium;  there  are  three  new  features  in  this  section.  First  the 
proof  of  Stoppelli's  results  is  considerably  simplified.  Second,  the 
results  are  global  relative  to  the  rotation  group.  Finally,  the  stabi¬ 
lity  of  the  solutions  is  determined.  The  number  of  solutions  will  be 
classified  by  load  types;  this  classification  scheme  is  explained  in 
Section  6.  (Some  work  related  to  the  "type  classification"  was  given  by 
Ogden  [1977]).  In  Section  7  a  second  order  bifurcation  equation  is  shown 
to  be  a  gradient.  This  consequence  of  Betti  reciprocity  is  basic  to  our 
analysis.  Section  8  gives  a  complete  bifurcation  analysis  of  type  1 
loads  (the  case  considered  by  Stoppelli),  including  a  stability  analysis. 
Mew  local  and  global  solutions  are  found.  The  final  section  makes 
explicit  the  comparison  with  Stoppelli's  theorem. 

The  second  paper  of  this  series  will  analyte  the  remaining  types 
2,  3  and  4,  using  a  reformulation  of  our  gradient  results,  discuss 
linearisation  stability,  parallel  loads  and  will  give  additional 

connections  with  the  literature. 
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Statement  of  the  Problem 


3  * 

Let  8  C  ir  be  an  open  bounded  set  with  smooth  boundary  and 

assume  for  convenience  that  0  £  8.  Let  1  <  p  <  *>,  s  >  (3/p)  +  1 

and  let  C  be  the  space  of  maps  .-p :  8  ->■  !R  that  are  of  class  W 

1  so 

(so  they  are  C  )  such  that  6(0)  =  0,  a  is  a  IV  ’^-diffeomor  '  sm 

onto  its  image,  and  J (ip )  >  0,  where  J(.j)  is  the  Jacobian  of 

For  example,  if  y:  8  ^  IRJ  is  W^’^-close  to  the  identity  and 

p(0)  =  0,  then  i  t  C  .  If  Q  is  a  linear  isomorphism  of  IR'3  to 

R-3  with  det  Q  >  0,  then  Qo>p  €  C  as  well. 

Let  points  in  8  be  denoted  X  £  8  and  points  in  RJ  be 

denoted  x.  Sometimes  we  write  x  =  *  ( X ) .  Let  T^.8  be  the  tangent 

space  to  8  at  X,  regarded  as  vectors  in  IR-3  based  at  X.  We  do 

not  identify  T^S  and  IR'3  for  conceptual  clarity.  For  j  £  C  ,  let 

F(X)  £  L(T^8,  JR3)  be  the  derivative  of  i  at  X;  by  standard  abuse 

of  notation  we  write  F(X)  =  Qp(X)  or  76 (X)  interchangeably  and 


L(T^8,  IR"3)  denotes  the  set 

of  all 

linear  maps 

of 

Txs 

to  IR'3 . 

We  let  F(X)T  £  L0RJ,  TXS) 

denote 

the  adjoint 

of 

F  (X) 

relative  to 

the  Euclidean  inner  product.  Observe  that  F(X)  £  L+(T^.8,  IRJ)  ,  the 
linear  transformations  with  positive  determinant,  since 
det  F (X)  =  J(p)(X)  >  0.  We  let  C  =  FTF  (that  is, 

CfX)  =  F(X)  F(X)  £  L(T^B,  )  denote  the  Cauchy-Green  tensor. 


We  believe  that  our  results  also  hold  when  8  has  piecewise  smooth 
boundary.  This  program  depends  on  elliptic  regularity  for  such  regions. 
Except  in  special  cases,  this  theory  is  non-existent  and  seems  to 
depend  on  a  modification  of  the  usual  Sobolev  spaces  near  corners. 
However,  for  simple  shapes  like  cubes,  the  necessary  regularity  can 
be  checked  by  hand  in  situations  where  the  linearised  elastostatic 
equations  can  be  solved  explicitly. 


4 


Observe  that  C(X)  £  Svm  (Tv8,  Tv8')  ,  the  positive  definite 

•  pos  XX  r 

symmetric  linear  transformations  on  T^B  . 

Assume  we  are  given  a  smooth  stored  energy  function  W  defined  on 

Dairs  CX. Cl  where  C  €  svm  (T„B,  T„S).  For  )£C  .  the  stored  energy 

•  pos  X  X  - t- 

of  j  is  !  W(X,C(X) )dV(X) ,  where  C  is  the  Cauehv-Green  tensor  of  3 
~  8 

and  dV  is  the  volume  element  in  S  .  The  fact  that  V  depends  on  3 

only  through  the  ooint  values  of  C  is  usually  called  material  frame 

indifference.  (See  Truesdell  and  Noll  [ 1965]  and  Marsden  and  Hughes 

[1978]  for  discussions.)  Since  C  is  a  function  of  F,  we  shall, 

T 

by  abuse  of  notation  also  write  W(X,F)  for  W(X,F  F) . 

The  first  Piola-Kirchhof f  stress  tensor  P(X,F)  is  defined 

3 

by  P(X,F)  =  ^-=-W(X,F),  the  partial  derivative  of  W  with  respect 

O  r 

to  F.  Thus,  P(X,F)  £  LCT^BJR'5)  .  The  second  Piola-Kirchhoff  stress 
tensor  S(X,C)  is  defined  by  S(X,C)  =  zlrHVfX.C),  so 

-  cf  L 

_  * 

S(X,C)  £  symCT^B,  T^B)  .  From  the  chain  rule,  one  has  the  relationship 
P(X, F) •  G  =  4s(X,C)-'FTG  *  gV 
for  all  G  £  LfT^BjR-5)  .  ' 

For  finite  dimensional  vector  spaces  V,W,  the  bilinear  map 
L(IV,  V)  x  L(V,IV)  -*•  IR;  (A,B)  h-  itrace(A  B)  defines  an  isomorphism 

■k 

L(V,IV)  =  L(1V,V).  If  V  =  IV  is  an  inner  product  space,  this 

* 

isomorphism  identifies  sym(V,V)  with  syn(Y,V).  Using  this  identifi¬ 
cation,  we  get  PfX,F)  £  L(1R'’,  T^8)  ,  S(X,C)  £  svmlT^BjT^B)  and 
P(X,F)  =  S(X,C)°FT,  or  P  =  SFT  for  short, 

“The  mass  density  does  not  appear  in  our  formulas  as  we  are  building 
it  into  the  definitions  and  work,  for  example,  with  hodv  force 
per  unit  volume  rather  than  per  unit  mass. 


o 


Let  A(X,F)  =  -<-=-l.X,F)  denote  the  elasticity  tensor.  Identifying 

a  r  — . —  —  -  -  ■  - 

L(LlT^S JRJ)  ,  l(T^S,  IR3)*)  in  the  usual  way  for  second  derivatives 
with  B(l(T^S,  IR3)),  the  space  of  bilinear  maps  of 
L(T^5,  IR3)  *  L(T^8,  IR3)  to  IR  ,  A(X,F)  determines  a  symmetric 
bilinear  map  on  L(T^B,  IR3)  . 

The  second  elasticity  tensor  C(X,C)  is  similarly  defined  to 


be  —  * 
e  3C  ~ 


3“W 


-rg evaluated  at  (X,C)  ,  so  may  be  regarded  as  a 
symmetric  bilinear  map  on  sym(T^3,  T^B) .  The  chain  rule  gives 


2A(X,F)  *(G,H)  =  C(X,C)-(FTH  *  HTF,  FTG  +  GTF)  +  S(X  ,C)  *  [HTG  +  GTH) 


The  following  two  assumptions  will  be  made  in  the  first  two 
papers  of  this  series. 

(HI)  The  undeformed  state  is  stress  free;  i.e.  P(X,I)  =  0, 
or  equivalently,  S(X,I)  =  0,  where  I  is  the  identity. 

(H2)  Strong  ellipticity  holds:  there  is  an  e  >  0  such  that 

A(.C,I)  *(v  ®  £,  v  &  £)  >  e !l C ii ~  11  v It - 

for  all  £  €  T^B  and  v  €IR3,  where  v  0  £  £  L(T  8,  IR3)  is  defined  hv 

'  (v  0  £)  (V)-v£  (\  ' . 

The  classical  elasticity  tensor  is  defined  by  C(X)  =  -C(X,I) , 
so  c(X)  is  a  symmetric  bilinear  mapping  on  sym(T^B,T^B)  to  IR ;  at 
<j>  =  I,  we  identify  T^B  and  IR"5  since  x  and  X  coincide.  By 


A(X,I) • (G,H)  =  jC(X)-(G  +  GT,  H  +  HT) 


(HI)  , 


6 


I 

! 

i  Regarding  A(X,I)  LCLiT^R.T^S) ,  L(T^8,  Tx8))  and 

C(X)  S  L(sym(TxB,  Tx«) ,  sym(Tx8,  T^B) )  .  this  reads 

2A(X,I)*G  =  d(X)  •  (G  ♦  GT)  , 

or,  if  G  is  symmetric, 

A(X,I)*G  =  c(X)(G)  • 

By  (H2) ,  solvability  of  the  linearized  equations  of  elastostatics  can 
be  determined  by  the  Fredholm  alternative  (see,  e.g.,  Marsden  and 
Hughes  [1978]). 

We  shall  let  B:  8  IR'5  denote  a  given  body  force  (per  unit 
volume)  and  T :  38  -  JRJ  a  given  surface  traction  (per  unit  area). 

These  are  dead  loads;  in  other  words,  the  equilibrium  equations 
for  that  we  are  studying  are: 


(E) 


DIV  P(X,  F(X))  +  B(X)  =0  for  X£  8 

P(X,F (X) ) -N(X)  =  T(X)  for  X  t  36 


where  N(X)  is  the  outward  unit  normal  to  38  at  X  €  36  and  DIV  F 
is  the  divergence*  of  P(X,F(X))  with  respect  to  X. 

Let  L  denote  the  space  of  all  pairs  (B,t)  =  l  of  loads  (of 
class  IV'  “’P  on  5  and  VIs  ^  on  3R  )  such  that 


j  B(X)dV(X)  +  t ( X ) dA ( X )  =  0 

‘Recall  from  above  that  P(X,F(X))  C  L(1R'\  Tx6) 
then  P(X,F(X))v  defines  a  vector  field  on  8  ; 
defines  the  vector  field  DIV  P  by  (DIV  P)-v  = 


For  any  v  fc  R-5  , 
its  divergence 
fHV  (Pv)  . 


i.e.  the  total  force  on  8  vanishes,  where  dV  and  dA  are  the 


respective  volume  and  area  elements  on  8  and  38.  Using  the 
divergence  theorem ,  observe  that  if  the  pair  (B,t)  is  such  that  (E) 
holds  for  some  )  6  C,  then  (B,t)  S  L  . 

Throughout  the  paper,  the  group  S0(3)  =  {Q  6  L(1R5JR3)[  Q^Q  =  I 
and  det  Q  =  +1}  of  proper  orthogonal  transformations  will  play  a 
key  role. 

By  (HI)  ,  =  I g  (the  identity  map  on  8)  solves  (E)  with 

B  =  t  =  0.  By  material  frame  indifference,  0  =  Q j  8  is  also  a 
solution  for  any  Q  €  SO(3) .  The  map  Q  h-  Q | B  embeds  SO (3)  into 
C  and  we  shall  identify  its  image  with  SO(3).  Thus,  the  "trivial" 
solutions  of  (E)  are  elements  of  SO(3) . 

Our  basic  problem  is  as  follows: 

(PI)  Describe  the  set  of  all  solutions  of  (E)  near  the  trivial 
solutions  SO(3)  for  various  loads  J  6  L  near  cero. 

Here  "describe"  includes  the  following  objectives: 

(a)  counting  the  solutions 

(b)  determining  the  stability  of  the  solutions 

(c)  showing  that  the  results  are  insensitive  to  small 
perturbations  of  the  stored  energy  function  and  the  load. 


This  section  is  devoted  to  the  geometry  of  the  load  space  L. 
Many  of  the  results  of  this  section  are  available  in  the  literature, 
but  we  gather  them  here  for  convenience. 

Before  beginning,  we  shall  recall  a  few  notations  and  facts 
about  the  rotation  group  S0(3) .  Let 

M,  =  L(R'>,  IR^)  =  linear  transformations  of  IR3  to  IR"5 

sym  =  {A  S  M.|AT  =  A} 

skew  =  {A  £  M,|AT  =  -A} 
j  1 

We  identify  skew  with  so(3),  the  Lie  algebra  of  S0(3);  skew  and 
are  isomorphic  by  the  mapping  v  SIR"5  n-  v  €  skew  ,  where 
v(w)  =  w  x  v.  If  v  =  (p,q,r)  relative  to  the  standard  basis,  then 


The  Lie  bracket  is  [v,w]  =  v&  w-  w®v  =  -(v  x  w)  where 
v  ®  w  €  M,  is  given  by  (v  ®  w) (u)  =  v<w,u>  .  The  inner  product 
is  <v,w)  =  itrace  (v  w) ,  the  Killing  form  on  so(3).  Finally, 
exp(v)  is  the  rotation  about  the  vector  v  in  the  positive  sense, 
through  the  angle  tl v II . 

Now  we  turn  to  a  study  of  L  .  For  f  €C  and  i€  L  ,  we  say 


that  i  is  equilibrated  relative  to  ?  if  the  total  torque  in  the 
configuration  f  vanishes: 


9 


|  4> (X)  x  BlX)dV(X)  +  f  <MX)  x  T(X)dAI.X)  =  0 
J5  '  38 

where  l  =  (B,x).  From  symmetry  of  the  stress  one  sees  that  if 
l  =  (B,t)  S  L  satisfies  (E)  for  some  t>  *=  C,  then  l  is  equilibrated 
relative  to  $.  (.An  easy  proof  uses  the  Piola  transform;  of.  Marsden 
and  Hughes  [1978]). 

Let  denote  the  loads  that  are  equilibrated  relative  to  the 

identity  configuration  Ig  . 

Define  the  astatic  load  map  k:  L  *  C  -*■  M.  by 

k (1,0)  =  B(X)  X><>(X)dV(X)  +  [  t(X)  ».i>(X)dA(X) 

J3  J38 

and  write  k(£)  =  k(£,  Ig) . 

We  have  actions  of  S0(3)  on  L  and  C  given  by: 


Action  of  S0(3)  on  L:  Q£(X)  =  (QB(X)  ,  Qt(X)) 
Action  of  S0(3)  on  C :  Q<J  = 


Note  that  Q9.  means  "the  load  arrows  are  rotated,  keeping 

the  body  fixed."  We  shall  write  0 ^  and  for  the  SO ( 3 )  orbits 

of  1  and  <fc.  Thus,  0  denotes  the  trivial  solutions  corresponding 

TB 

to  1  =  0. 


The  following  is  3  list  of  basic  observations  about  the  astatic 


load,  each  of  which  is  readily  verified: 


10 

(Al)  l  is  equilibrated  relative  to  f  if  and  only  if 

k(£,$)  G  sym.  In  particular,  l  G  Lg  if  and  only  if 
k (l)  G  sym. 

(A2)  (equivariance) .  For  l  G  L ,  )  G  C.and  ,Q_,  G  S0(3) 

k(Qii,Q2.5)  =  Q1k(i,i)Q]1 
In  particular,  k(Q£)  =  Qk(2) 

(A3)  (infinitesimal  equivariance).  For  £GL,  d>  G  C  ,  and 
W  W.,  G  skew, 

k(WjJl^)  =  WjkU.O),  k(^ ,W2d>)  =  -k(il,?)  W, 

In  particular,  k(WJL)  =  Wk (2) . 

Later  on,  we  shall  be  concerned  with  how  the  orbit  of  i  G  L 
meets  L  .  The  most  basic  result  in  this  direction  is  the  following. 

3.1.  DaSilva's  Theorem.  Let  l  G  i.  Then  0„  -'A  L  /  Id. 

Proof.  By  the  polar  decomposition,  we  can  write  k(2) 

for  some  Q  G  S0(3)  and  A  G  sym.  By  LA2)  ,  k(Q2)  =  QkiL’i 

so  by  (Al)  ,  Q£  G  L  .  ■ 
e 

Similarly,  any  load  can  be  equilibrated  relative  to  any  chosen 
configuration  by  a  suitable  rotation. 

The  concept  of  an  axis  of  equilibrium  deals  with  the  case  in 
which  C  meets  L,  in  a  degenerate  way. 


T 

Q  A 


=  A  G  svm, 
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3.2.  Definition.  Let  i  €  L  and  v  6  1R'5,  llvll  =  1.  We  sav  that 

-  e 

v  is  an  axis  of  equilibrium  for  l  when  exp(9vK  e  L  for  all 
real  9,  i.e.  when  rotations  of  l  about  the  axis  v  do  not 
destroy  equilibration  relative  to  the  identity. 

A  number  of  useful  ways  of  reformulating  the  condition  that  v 
be  an  axis  of  equilibrium  are  as  follows. 

3.5.  Proposition.  Let  i  €•  L and  A  =  k(2)  €  sym.  The  following 
conditions  are  equivalent: 

1 .  Z  has  an  axis  of  equilibrium  v 

2.  there  is  a  v  E  1RJ  ,  llvll  =  1  such  that  v£  €  Lg 

3.  W  h-  AW  +  WA  fails  to  be  an  isomorphism  of  skew  to  itself 

4.  trace  A  is  an  eigenvalue  of  A. 

Proof.  1  =»  2.  This  follows  by  differentiating  expi~v).'.  in  9 
at  9  =  0.  2  =>  1.  By  fA2)  , 

k(exp(9v)J,)  =  ‘I  +  (9v)  +  -r’lOv)-  +  .  ..'k(il)  . 

Since  k(v’)  =  vk('l)  is  symmetric,  this  is  symmetric,  term 
by  term. 

2  =*  3.  Since  k(vL)  =  v.A  is  symmetric,  vA  +  Av  =  0,  so 
W  K  AW  +  WA  is  not  an  isomorphism. 

3^2.  There  exists  a  v  £  ir\  llvll  =  1,  such  that  vA  +  Av  =  0, 
so  kfvL)  =  vA  is  symmetric. 

3  «*  4 .  Define  L  £  M.  bv  L  =  (trace  AjI  -  \.  Then  one  has 

.> 


the  relationship 
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(Lvr  =  Av  +  vA 

(In  fact,  if  [u,v,w]  denotes  the  triple  product,  the  relationship 

[Bu,Bv,Bw]  =  det  B[u,v,w]  gives  [Au.v.w]  +  [u,Av,  v;]  +  [u,v,Aw] 

T 

=  (trace  A)  [u,v,wj  .  This  yields  (Lv)  =  A\A  +  A  v  ,  which 
gives  the  claimed  results  for  symmetric  A).  Therefore, 

Av  +  vA  =  0  if  and  only  if  Lv  =  0  ,  i.e.,  v  is  an  eigenvector 
of  A  with  eigenvalue  trace  A  .  • 

3.4.  Corollary.  Let  l  £  and  A  =  k (£)  G  svm.  Let  the  eigen¬ 
values  of  A  be  denoted  a,  b,  c.  Then  l  has  no  axis  of  equili- 
brium  if  and  only  if  (a  +  b) (a  +  c) (b  +  c)  /  0. 

Proof.  This  condition  is  equivalent  to  saying  that  trace  A 
is  not  an  eigenvalue  of  A.  B 

3.5.  Definition.  We  shall  say  that  i  G  Lg  is  a  type  0  load  if  l 
has  no  axis  of  equilibrium  and  if  the  eigenvalues  of  A  =  k(£) 
are  distinct. 

The  following  shows  how  the  orbits  of  type  0  loads  meet  Lg. 

5.6.  Proposition .  Let  i  £  Lg  be  a  type  0  load.  Then  Q  Cl  L? 
consists  of  four  type  0  loads. 

Proof .  We  first  prove  that  the  SO(5)-orbit  of  A  in  M_  under 
the  action  0-QA  meets  sym  in  four  points.  Relative  to  its  basis 
of  eigenvectors,  we  can  write  A  =  diag(a,b,c).  Then  0 ^ 
contains  the  four  points 


diag(a,b,c) 
diag(-a,-b,cj 
diug(-a,b , -c) 
diag(a,-b,-c) 


(Q  =  I) 

(Q  =  diagf-1,-1,1)) 
(Q  =  diagC-1,1-1)) 
(.Q  =  diag(l,-l , - 1 ) ) 


These  are  distinct  points  since  (a  +  b) (a  +  c)  (b  +  c)  i  0.  Now 
suppose  a,  b  and  c  are  distinct.  Suppose  QA  =  S  €  sym.  Then 

S"  =  A-.  Let  be  an  eigenvalue  of  S  with  eigenvector  u. . 

2  g  a  ■> 

Then  S  u^  =  y7u^  =  ^~u- >  so  is  an  eigenvalue  of  A".  Thus, 

2 

as  the  eigenvectors  of  A“  with  a  given  eigenvalue  are  unique, 
u^  is  an  eigenvector  of  A  and  -y  is  the  corresponding  eigenvalue. 
Since  det  Q  =  +1,  det  S  =  det  A,  so  we  must  have  one  of  the  four 
cases  above. 

By  equivariance ,  k(O^)  n  sym  =  ft  sym  is  a  set  consisting  of 

four  points.  Now  0 ^  T,  Lq  =  k- ^  ~  ’  so  suffices  to 

show  that  k  is  one-to-one  on  0^.  This  is  a  consequence  of  the 
following  and  (A2) . 


3.7.  Lemma .  Suppose  A  €  sym  and  dim  ker  A  <  1.  Then  A  has  no 
isotropy :  i , e .  QA  =  A  implies  Q  =  I . 


Proof.  Every  Q  /  I  acts  on  1R'  by  rotation  through  an  angle 
0  about  a  unique  axis,  that  is,  about  a  line  through  the 
origin  in  IR’’  .  Now  QA  =  A  means  that  Q  is  the  identity  on  the 
range  of  A.  Therefore  if  0^1  and  QA  =  A,  the  range  of  A 
must  be  0  or  1  dimensional,  i.e.  dim  ker  A  J.  ■ 
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Finally  in  this  section,  we  study  the  range  and  kernel  of 
k:  L  -*  M- 

o 

3.8.  Proposition.  1.  ker  k  consists  of  those  loads  in  L for 
which  every  axis  is  an  axis  of  equilibrium. 

2.  k:  L  -*•  M-  is  surjective. 

Proof.  1.  Let  ?,  S  ker  k.  For  W  6  skew,  k(W)t)  =  Wk(£)  =  0 

so  WP.  £  Le',  by  3.3  every  axis  is  an  axis  of  equilibrium.  Conversely, 

if  W Si  €  L  for  all  W  €  skew,  then  k(WS.)  =  WkU)  is  symmetric 
e 

for  all  W;  i.e.  k(P,)W  +■  Wk(P)  =  0  for  all  W.  From 

(Lv)~  =  Av  +  vA,  where  A  =  k(P,)  ,  and  L  =  (trace  A) I  -  A,  we 

see  that  L  =  0.  This  implies  trace  A  =  0  and  hence  A  =  0. 

To  prove  2  introduce  the  SQ(3)  invariant  inner  product  on  L: 

t 

<l,i)  =  <  B(X),B(X))dV(X)  +  |  <T(X),T(X)>dALX). 

8  '  38 

T 

Relative  to  this  and  the  inner  product  trace (A  B)/2  on  M,,  the 
T 

adjoint  k  :  M  L  of  k  is  given  by 

kT I A j  =  (B,t)  where  BfXj  =  AX-G,  rfX)  =  AX-G, 

and 

G  =  f  f  AX  dVfX)  +  f  AX  JA  ( X 1  "^  /  !  j  GV  +  '  JA  . 

‘-'8  J';E  J  -(8  -38  J 

If  kT(A)  =  (0,01  then  it  is  clear  that  A  =  0.  It  follows  that 
k  is  surject i ve.  B 


I 
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3.9.  Corollary .  1.  ker  k  is  the  largest  subspace  of  Lg  that  is 

S0(3)  invariant . 

2.  k]  (ker  k)-*-  :  (ker  k)^  -*■  M,  is  an  isomorphism. 

Let  j  =  (k| (ker  k)^)'1  and  write 

Skew  =  j  (skew)  Sym  =  j  (sym) 


These  are  linear  subspaces  of  L  of  dimension  3  and  6  respectively. 
Thus  we  have  the  decomposition: 


i 


S0(3)  invariant  pieces- 


L  =  Skew  (to  Sym  ®  ker  k 


L 

e 


corresponding  to  the  decomposition  M.  =  skew  sym; 

1  T  1  T 

U  =  -  U  )  +  j(U  +  U  ) . 


.Mote:  Skew  and  L  need  not  be  orthogonal. 
-  e 


J4.  Equivalent  Reformulations  of  the  Problem 


Define  $:  C  -*  L  by  i(?)  =  (-DIV  P,  P - \T)  i.e. 


$($)(X)  =  (-DIV  P(X,F(X)),  P (X , F CXI )  *NI (X) ) 


so  the  equilibrium  equations  (E)  become  f($)  =  L  .  From  material 
frame  indifference  we  have  equi variance  of  i;  i.e.  HQc)  =  QMi), 
Standard  Sobolev  estimates  show  that  i  is  a  smooth  mapping  (see, 
for  example,  Palais  [1968]).  The  derivative  of  5  is  given  by 


D<f((J))*u  =  (-DIV (A« Vu) ,  (A*7u)  »N) 


and  at  $  =  IR  this  becomes 


D$(I  )«u  =  (-DIV(c*e),(c*e)»N) 


1  T 

where  e  =  —  [7u  +  (Vu) 

I f  D4>(Ir)  :  T.  C  ■*  L  were  an  isomorphism  we  could  solve 

8  8 

$($)  =  i  uniquely  for  i>  near  Ig  and  l  small.  The  essence 

,  .  ,  .  rsince  ?CS0(3)  )=0 ,  ■ 

of  our  problem  is  that  D*(Ig)  is  notion  isomorphism,  -  D^TT'contair 

Define  C  =  fu  £  Tt  C1  u(0)  =  0  and  7u(0)  £  svm}  .  From 
sym  • 

o 

(H2)  and  linear  elasticity,  we  have: 


4.1.  Lemma.  CK*>  T I  _)  I C  :  C  -1-  L  is  an  isomorphism. 
-  8  1  svtti  svm  e  - 1 - 


I 


The  connection  between  the  astatic  load  map  k :  L  —  M_  and  i 
is  seen  from  the  following  computation  of  k-$. 

-I  2 .  Lemma .  Let  v  €  C  and  let  P  be  the  first  Piola-Kirchnoff 
stress  tensor  of  5  .  Then 

k  O i  : ) )  =  ’  pdV 

■;S 

This  follows  by  an  application  of  Gauss'  theorem  to 

r  r 

=  |  f-DIV  P)  »  X  dV (XI  +  '  (P  N;  S>  X  dAfXl. 

J  6  j  38 

This  should  be  compared  with  the  astatic  load  relative  to  the 
configuration  J>  rather  than  Ig;  one  gets 

k($(2>) ,?)  =  !  S  dV 

J8 

which  is  symmetric,  while  k(1>(£))  =  k(l5(b),I_)  need  not  be. 

o 

To  study  solutions  of  $($)  =  i  for  <J>  near  the  trivial 

solutions  and  l  near  a  given  load  2^,  it  suffices  to  take 

.  This  follows  from  Da  Silva's  theorem  and  equivariance  of 

Let  (-svm  be  regarded  as  an  affine  subspace  of  C  centered  at 

I„.  Let  ?  be  the  restriction  of  T  to  C  .  From  the  implicit 
8  sym  r 


function  theorem  we  get: 
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Later  we  shall  show  how  to  compute  D~F(0)  in  terms  of  ( I g) 
and  C. 

Now  we  are  ready  to  reformulate  problem  (PI). 

(P2)  For  a  given  G  near  zero,  study  how  0 ^  meets  the  graph 
of  F  for  various  2.  near  . 

Problems  (PI)  and  (P2)  are  related  as  follows.  Let  o  solve  (E) 
with  i  G  L  and  Q  be  such  that  $  =  Q<J>  £  Then  $(5)  =  Q2.  , 

so  the  orbit  of  l  meets  the  graph  of  F  at  $($).  Conversely,  if 
the  orbit  of  j l  meets  W  at  4>(<j>)  =  Q£  ,  then  <+>  =  Q  ^ b  solves  (E) . 
We  claim  that  near  the  trivial  solutions,  the  numbers  of  solutions 
to  each  problem  also  correspond.  This  follows  from  the  next  lemma. 


4.4.  Lemma.  There  is  a  neighborhood  U  of  Ig  in  C  such  that 
0  G  U  and  QcJ)  G  U  implies  Q  =  I  g. 


Proof.  Note  that  C  is  transverse  to  C  at  I R  and  IK 
-  sym  I  g  0 

has  trivial  isotropy.  Thus,  as  S0(5)  is  compact,  0  is  closed,  so 

S 

there  is  a  neighborhood  of  Uq  of  Ig  in  C_vm  such  that 

Q| B  G  UQ  implies  Q  =  I.  The  same  thing  is  true  of  orbits  passing 


through  a  small  neighborhood  of 


I ^  by  openness  of  transversalitv  and 


compactness  of  S0(3) .  ■ 


If  0  meets  M  in  k  points  =  +>(uk),  i  =  l,...,k  then 

S.  are  distinct  as  4>  is  1-1  on  a  neighborhood  of  I„  in  C 
i  b  sym 

If  this  neighborhood  is  also  contained  in  U  of  4.4,  then  the  points 


:o 


Qi  4k  =  2k  are  also  distinct  by  l.l.  Hence  the  problems  l  PI')  and  ;  P  J 
are  equivalent. 

In  connection  with  the  action  (Q,A)  h-  ')A  of  S0(3)  on  M.  we 
shall  require  some  more  notation.  Let 

Skew  (A)  =  -  A^l  s  skew 

and 

Sym(A)  =  i(A  +  A^l  ^  sym 

be  the  skew  symmetric  and  symmetric  parts  of  A,  respectively. 

We  shall,  by  abuse  of  notation,  suppress  j  and  identify  Sym 
with  sym  and  Skew  with  skew.  Thus  we  will  write  a  load  l  £  L  as 
i  =  (A,n)  where  A  =  k(£)  t  M.  and  n  6  ker  k;  hence  1  6  1^  precisely 
when  A  €  sym.  The  action  of  SO  (3)  on  L  is  given  by 

Qe  =  (QA, Qn). 

Using  this  notation  we  can  reformulate  problem  (P2)  as  follows: 

(P3)  For  a  given  l ^  =  (A^.n^)  €  L e  near  cero,  and  l  =  (A,n) 
near  £  find  Q  €  SO (31  such  that 

Skew(QA  F(Sym(QAl ,Qnl  =  0. 

Define  the  rescaled  map  F :  IR  <  L  -*■  Skew  by 


(3.2a) 


(3.2b) 


:i 


F(A,2)  = 

A~ 

Since  F (.O')  =  0  and  DF(O)  =  l),  F  is  smooth.  Moreover,  if 
F(Z)  =  4g(‘0  +  ^C(i)  +  ...  is  the  Taylor  expansion  of  F  about 
zero,  then  F ( \ ,i )  =  +  ^-C ( l)  +  ...  . 

In  problem  (F.)  let  us  measure  the  size  of  L  by  the  parameter 
A  .  Thus,  replace  =  i  for  1  near  zero  by  >($)  =  \l 

for  A  near  zero.  This  scaling  enables  us  to  conveniently  distinguish 
the  size  of  i  from  its  'orientation'.  In  the  literature  l  has 
always  been  fixed  and  A  taken  small.  Here  we  allow  L  to  vary  as 
well.  Thus  we  arrive  at  the  final  formulation  of  the  problem. 

(P4)  For  a  given  =  (Afl,n^)  <=  ,  for  l  near  and  A  small , 

find  Q  €  S0(5)  such  that 

Skew(QA)  -  AF(\,  Sym(QA),Qn)  =  0. 

The  left  hand  side  of  this  equation  will  be  denoted  H(A,A,n;Q) 


or  H(A,Q)  if  A,  n  are  fixed. 


We  shall  begin  the  analysis  by  giving  an  (almost  trivial)  proof 


of  one  of  the  basic  theorems  of  Stoppelli  ’ 1938) 

5.1.  Theorem.  Suppose  l  <=  has  no  axis  of  equilibrium.  Then  for 
a  sufficiently  small,  there  is  a  unique  o  £  ^svm  and  a  unique  Q 
in  a  neighborhood  of  the  identity  in  S0(3)  such  that  p  =  Q~*i>  solves 
the  traction  problem 

H<p)  =  U  . 

Proof.  Define  H:  1R  x  S0(3)  -*•  Skew  by 

H(X ,Q)  =  Skew(QA)  -  XF (X ,Sym(QA) .  Qn) ) 

where  i  =  (A,n)  £  =  Sym  ker  k  is  fixed.  Mote  that 

D,H(0 , 1)  •  W  =  Skew  (  WAJ  =  *  AW).  By  Proposition  3.3,  this  is 

an  isomorphism.  Hence,  by  the  implicit  function  theorem,  H(X,Q)  =  0 
can  be  uniquely  solved  for  Q  near  I  £  S0(3)  as  a  function  of  \ 
near  0  €  IR.  ■ 

The  geometric  reason  "why"  this  proof  works  and  the  clue  to 
treating  other  cases  is  the  following. 


The  only  other  com;)  lete  proof  in  F.nglish  we  know  of  is  given  in 
Van  Buren  ;1968),  although  sketches  are  available  in  (Irioli  19t>3  .  , 
Truesdell  and  Noll  )19bS)and  Wang  and  Truesdell  19”3)  .  Our  proof 
rather  different;  the  use  of  the  map  P  avoids  a  series  of 
licated  estimates  used  by  Stoppelli  and  Van  Buren. 


is 
comp 


3  .  - 


Lemma.  A  load 


l  6  has  no  axis  of  equilibrium  precisely 

when  i.  =  *  T^O.,  .  In  particular,  if  l  has  no  axis  of  eauili- 

brium,  then  C,  intersects  L  transverselv  at  l 
-  I - e  - ‘ - 

Proof.  The  tangent  space  to  0^  at  16  Lg  is  T,d„  =  .  WJl 'W  6  skew; 
and  the  projection  of  this  into  the  complement  Skew  to  Lg  is 
\U  ’ > (WA  *  AW)  where  A  =  k(i)  .  The  result  then  follows  from 
part  3  of  3.3.  ■ 


We  have  shown  that  there  is  only  one  solution  to  £>(?)  =  M  near 
the  identity  if  A  is  small  and  l  has  no  axis  of  equilibrium.  How 
many  solutions  are  there  near  the  trivial  solutions  S0(5)?  As  we 
shall  see,  this  problem  has  a  non-trivial  answer  which  depends  on  the 
type  of  1.  We  analyte  the  simplest  case  here.  Recall  from 
definition  3.3  that  a  load  1  £  L ^  is  said  to  be  of  type  0  if  1  h3S 
no  axis  of  equilibrium  and  if  A  =  k(ll  has  distinct  eigenvalues. 

Loads  with  no  axis  of  equilibrium  occur  amongst  other  types  of 
loads  classified  in  the  next  section,  and  Stoppelli's  theorem  5.1 
applies  to  them.  However,  the  global  structure  of  the  solutions 
("global"  being  relative  to  S0(3))  is  different  for  the  different  types. 
For  type  0  the  situation  is  as  follows. 

5.5.  Theorem.  Let  €  Le  be  of  type  0.  Then  for  «  sufficient ly 
small,  f(t>)  =  \Jn  has  exact lv  four  solutions  and  •  in  a 

- _  '  0  - : -  1  a - ■! - 

neighborhood  of  the  trivial  solutions  SO ( 3 )  C  C  (see  Figure  5'. 


Proof.  By  3.6,  0,  ,  meets  in  four  points.  By  5.1,  in 

^0  6 

a  neighborhood  of  0  in  L,  0-.  .  meets  ^  in  exactly  four  points, 
the  images  of  5  ,  b,  and  5^,  say.  Thus  problem  (P2)  has  four 

solutions.  By  the  equivalence  of  (PI)  and  (P2) ,  so  does  (PI).  ® 


50(3) ;  the 
trivial  solutions 


Let  A  =  k[£.g)  and  =  { Q I QA  G  svm } .  From  tlie  proof  of  5 .  t> 
we  see  that  5^  is  a  four  element  subgroup  of  SO  (S')  isomorphic  to 
2-,  ®  Zr  By  our  earlier  discussions,  the  elements  5  are  obtained 
from  b ^  by  applying  rotations  close  to  elements  of  S^.  In 
particular,  as  \  -*•  0,  the  solutions  {p.|  converge  to  the  four 
element  set  (regarded  as  a  subset  of  C). 

For  i  sufficiently  close  to  ,  the  problem  ?(•:■)  =  will 
also  have  four  solutions.  Indeed  by  the  openness  of  t ransversa 1 ity , 

Ci,£  will  also  meet  SI  in  four  points.  In  other  words,  the  picture 
for  type  0  in  Figure  2  is  stable  under  small  perturbations  of 
Next  we  study  the  stability  of  the  four  solutions  found  in 
Theorem  S.5.  This  will  be  done  under  the  hypothesis  that  the  classical 
elasticity  tensor  in  stable.  IVe  introduce  the  following  condition: 

(Ho)  Assume  there  is  an  n  >  0  such  that  for  all  e  £  SymlT^S,  T^8) , 

e(e)  =  -c(X){e,e)  >  n*lell“,  II 'il  =  pointwise  norm 

( c(e )  is  the  stored  energy  function  for  linearized  elasticity). 

Because  of  difficulties  with  potential  wells  and  dynamical 
stability  in  elasticity  (see  Knops  and  Wilkes  (19"2]  and  Ball,  Knops 
and  Marsden  [ 19T8]  )  we  shall  adopt  the  following  "energy  criterion" 
definition  of  stability. 

5.4.  Definition .  A  solution  0  of  will  be  called  stable 


if  t> 


is  a  local  minimum  in  C  of  the  potential  function 


2b 


V,  i  ;)  =  W i  d ) dV  -  (l,;) 

) 

S 

r  r 

where  B(  X) -d  f.XIdV  (X)  +  t  f  X)  •  ;  (X)  dA(  X)  =  trace  ki'i,*,: 
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If  o  is  not  stable,  its  index  is  the  dimension  of  the  largest 
subspace  of  vectors  u  tangent  to  C  at  D  with  the  property 
that  ;  decreases  along  some  curve  tangent  to  a.  (Thus,  index  0 
corresponds  to  stability.) 


5.5.  Theorem .  Assume  (HI) -0151  and  let  i. be  as  in  5.5.  For 
sufficiently  small,  amongst  the  four  solutions  given 

— . - - - - -  1  -  O  4  - 

bv  5.5,  exactly  one  is  stable;  the  others  have  indices  1,  2,  and  3. 
Suppose  i  is  a  solution  approaching  Q  €  a£  \  -  0.  Then  * 
is  stable  if  and  only  if  QA  -  trace (QA) I  €  sym  is  positive-definite. 
In  general,  the  index  of  }  is  the  number  of  negative  eigenvalues 
of  QA  -trace  (QA) I  . 


Proof.  bet  >(1  ■=  C  solve  Me)  =  '• .  (,s  •' .  Then  is  a  critical 

point  of  V.  ,  .  Consider  the  orbit  c\  =  ;Qv(,  Q  ~  SO(5i  of  ;  . 

''-0  :n 

Its  tangent  snace  decomposes  T,  C  as  follows: 

•Q 


T.  C  =  T,  i\  -f  IT.  0. 

• o  'o  *0  'o  *n 


first  consider  V  restricted  to  IT_  0,  1*.  I 

'  •(!  Ml 


at  r  in  the  direction  of  i  -  iT,  Cg  ■  ~  i * 

Ml  Ml 


o'cond  ler: .at ;  v- 


,;-f:  -i 


■  u , 


V  Q|s 


this  becomes 


[  c(X)  •  (e(X)  ,e(X)  jdV(.X')  , 

8 


where 


1  T 

e  =  — (7u  +  (7u)  ).  This  is  larger  than  a  positive  constant  times  the 
? 

L~  norm  of  e,  by  (H3)  .  However,  since  u  is  in  (T_t\l-, 

^  j 

2  2 

Hell'.,  >  (constant)  flu  11“ .  by  Korn’s  ineauality  (see  Fichera  ’1972]). 

L"  H 

By  continuity,  we  have 


D“V. 


'*0  0 


(u,u)  ^  5 Bull" 
H 


if  u  is  orthogonal  to  0  at  in  ami  \  is  small.  This  implies 

*0  0 

that  s  is  a  minimum  for  V, ,  in  directions  transverse  to  C 

0  UQ  P0 

(Actually  one  can  see  that  ■;>  is  a  local  minimum  in  the  topology 

I 

of  C  on  (T  0  )“  by  using  Tromba's  "29 “6 ’  version  of  the  Morse 

*0  50 

lemma. ) 

Next,  consider  V. ,  restricted  to  C  .  Bv  material  frame 

0  J0 

indifference,  W  is  constant  on  t\  and  so  as  ;■  must  be  a 

critical  point  for  V  restricted  to  04  ,  it  is  also  a  critical 

'0  "0 

point  for  a Zn  =  l  restricted  to  C  (where  <  1 1 1  =  It 

0  >0 

suffices  therefore  to  determine  the  index  of  '  t\  at  b^.  The 

*  0 

result  is  now  a  consequence  of  continuity  and  the  limiting  case 
■  -  0  given  in  the  following  lemma  about  type  0  loads. 


5.n.  Lemma .  let  l  he  type  0  and  let  A  =  k (  ?. )  .  Then  S  , 

regarded  as  a  subset  of  C  equals  the  set  of  critical  points  of 

..  C1  r  .  These  4  critical  points  are  nondegenerate  wit.n  indices 

8 

0,  I,  2  and  5;  the  index  of  Q  is  the  number  of  negative  eigenvalue 


of  t race 'OA) I . 


Proof.  First  note  that  L  =  (T.  SOlilT  since  D4>(I„1  has 
■  e  1^  o 

o 

kernel  T  S0f5)  =  skew,  lias  range  L  and  is  se lf-adi oint .  Thus 

*■  -o 

o 

Qi  £  L  if  and  onlv  if  l  i  T  S0t5) .  It  follows  that  Q£  €  L  if 
e  T  x  e 

j  0 

and  only  if  Q  is  a  critical  point  of  i  0  .  Recall  that 

1  8 

elements  of  =  {Q  c  St'f5):Q2  £  L  }  are  symmetric,  j 

To  compute  the  index  of  ,,  0  at  Q  e  S  we  compute  the 

8 

second  derivative 


t=0 


Now 


l(W"Q)  =  trace  k(l,lV'Q)  =  trace[kf ,ni»' 


.  i  -> 

=  trace  ;  AQ  w~  ]  =  trace  [WQA]  ' 


because  Q  =  Q  .  Tins  cmadratic  form  on  skew  is  represented  bv  the 
element  OA  -  trace (QA) 1  of  sym  as  is  seen  from  vA  +  Av  =  <  l.v  I  ’  with 
A  replaced  by  OA  and  trace  ivwl  =  dvw  .  "sing  the  representations 
for  fQA-  given  Proposition  3.t>,  namely 

Jiagla.b.ci,  di  ag  < -a  , -b  ,  c  '•  ,  di  ag( -a  ,h , -c  1  anti  di  ag  (  a  , -b  , -c  ) 


one  checks  that  all  four  indices  occur.  H 


Remark.  This  lemma  is  a  ~ncci.il  case  of  the  general  problem 


to  study  the  critical  points  of  near  *"unot  ional  <  on  orbits  of 


a 


i 


■’ll 


representation  of  a  Lie  group.  This  situation  will  arise  again  in 
our  analysis  of  the  other  load  types,  of.  Frankel  [1965]  a;u- 
Ramanujam  _19691. 
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it.  Class i fication  of  Orbits  in  M„ 

- o 

The  purpose  of  this  section  is  to  classify  orbits  in  M_  under 
the  action  (Q,A)  v*  qa  of  S0(3J  on  M_  by  the  way  the  orbits 
meet  svm.  It  is  enough  to  consider  orbits  0 of  elements  of  sym 
by  the  polar  decomposition.  We  begin  by  recalling  a  result  already 
proved. 

b.l.  Proposition  (Type  0).  Suppose  A  £  sym  has  no  axis  of  equilibrium 
and  has  distinct  eigenvalues.  Then  C1.  ~  sym  consists  of  four  points 
at  each  of  which  the  intersection  is  transversal. 

This  was  Proposition  3.6.  (Another  proof  of  this  is  given  below.) 

We  shall  let  the  eigenvalues  of  AS  sym  be  denoted  a,b,c. 

Using  the  terminology  from  53,  we  say  that  A  has  no  axis  of  equilibrium 
when  (a+b) (b+c) (a+c)  f  0;  i.e.  a  +  b  +  c  t  a,  b  or  c ,  and  in  this 
case  (?  intersects  sym  transversely  at  A. 

6.2.  Definition.  W'e  shall  say  A  is  of  type  1  if  A  has  no  axis 
of  equilibrium  and  if  exactly  two  of  a,b.c  are  equal  and  non-zero 

i say  a  =  b  ^  c ,  a^O). 

6.3.  Proposition .  I f  A  is  type  1,  then  C?  'A  sym  consists  of 

two  points  (each  with  no  axis  of  equi  librium!  and  an  Mr'  ^  ^  '  erica 

point  of  which  has  one  axis  of  equilibrium'1  . 

Before  proving  this,  we  give  a  number  of  lemmas  of  genera! 
utility.  If  l  t  1RP-  is  a  line  through  the  origin  in  IP  ,  let  0 
be  the  rotation  through  angle  ~  about 


I 


.  -I  .  Lemma  . 


Q,  is  an  embedding  of  5RP-  onto  501.5] 


Proof.  It  is  clear  that  I  h-  Q,  is  a  one-to-one  nap  of  FP" 

■>  ^  J  T* 

into  30(3).  Since  QT  =  I ,  Q,  =  Q,  =  Q,  so  it  maps  into 


SO  1.3)  sym. 


Every  Q=  SO <  3 )\ I  is  a  rotation  through  some  angle  5  about 


some  axis 


i.  If  such  Q  also  is  symmetric  then  it  has  three 


inderiendent  real  eigenvectors.  Hence  5 


0.5.  Coro  1 lary .  The  orbit  of  the  identity,  0 , .  meets  sym  m 
one  point  (I)  and  KP  =  (SO (3)  '  sym)\  I. 


6.6,  Lemma.  Let  A  €  sym  with  dim  kor  A  s  1  and  suppose  that 
Q  €  SO ( 3 V  I  and  QA  €  sym.  Then  0  =  0,  for  some  line  invariant 

under  A  ,  and  in  narticular  0  •=  sym. 


Proof.  We  can  suppose  Q  i  I.  There  is  a  unique  : up  to  si  erf 

unit  vector  x  E  IRJ  such  that  Qx  =  x.  Since  QA  -2  sym,  we  have 
7’ 

Q\  =  AQ  ,  so  QAQ  =  A.  Thus  QAx  =  Ax.  so  \x  =  ox  for  a 
constant  c.  Each  of  Q  and  \  leaves  \  =  x  ,  the  orthogonal 
complement  of  x,  invariant,  and  \  is  not  identicallv  zero  m  Y 

l.et  S  =  QA  t.  syn,  so  5“  =  A"  .  Since  Q’\’  is  a  rotation 

this  implies  S : V  =  Q.  A 1 V  giving  0=1  or  1  =  Q,  where 

i  f  x  l 

i.x'  is  the  line  through  x  .  Then  n  •=  Syn  as  in  Lemma  t» .  I  .  ■ 


4 


It  follows  t 
QA  sym,  then 
A  and  Q  can  be 


hat  if  dim  ken  \  1  im! 

QA  =  AQ,  so  as  0  is  both  orthogonal  and  symmetric, 
simultaneously  diagonal: ted. 


Proof  of  o.l.  If  A  has  u.j:inct  eigenvalues,  its  eigenvectors 
are  unique,  so  0  is  either  the  identity  or  is  a  rotation  by  ~  about 
one  of  the  eigenvectors.  ■ 

Proof  of  r.3.  Suppose  0  ?  a  =  b  *  c ,  and  let  w  be  an  eigenvector 
corresponding  to  the  eigenvalue  c.  Let  V  be  the  plane  orthogonal 
to  w  s  ;•  V  is  the  eigensp.ace  with  eigenvalue  a.  As  Q,  A  can  be 
-  mvoi.i  ■  ■■  diagonal! ted  and  Q  is  a  rotation  by  tt  (excluding 

'  =  .a-.e  either  Q  =  Q  or  0  *  Q,  for  l  a  line  in  V. 

•,  {  W  )  K, 

firmer  case,  Q.  ,  has  eigenvalues  (-a,-a,c)  so  has 
no  i\;s  of  euui  1  ibriuns.  In  the  latter  case ,  Q„A  has  eigenvalues 
a ,  -  a  ,  -  c  ;  ,  sc  w  is  an  axis  of  equilibrium,  (see  3.3).  B 

t>.~.  Corollary.  The  HT^  in  Proposition  6.3  is  a  right  coset  of 
the  subgroup  S*  of  all  rotations  about  w;  in  fact  HP ^  =  5_Q - 
where  ?,  is  any  fixed  line  in  V,  the  plane  orthogonal  to  w. 

Proof.  HP1  =  {Q?  |  ?,£  V}  and  we  have  the  easily  verified 


identity 


Q.,  =  exntvw  .Q. 

' i/2  '  "0 

where  i,  rakes  an  angle  -  with  in  the  Positive  sense  in  V  . 

.  -I 

These  lemmas  also  enable  us  to  handle  the  next  type. 

6..S.  Def  in  it  ion .  IVe  shall  say  A  is  of  type  1  if  \  has  no 

axis  of  equilibrium  and  all  three  of  a.b.c  are  equal  (and  so  4  0). 

6.9.  Proposition.  If  A  is  type  2,  then  C\  '  s>tti  consists  of 

"1 

one  point  (A)  and  an  ITOP-, 

Proof.  This  is  immediate  from  6.5.  9 

Notice  that  each  point  of  the  1RP~  has  a  whole  circle  of  axes 
of  equilibria;  namely  Q^A  has  as  axes  of  equilibria  all  vectors 
orthogonal  to  l.  The  eigenvalues  of  Q.,A  are  a, -a, -a. 

Types  0,1,  and  2  exhaust  all  symmetric  matrices  with  no  axis 
of  equilibrium,  and  it  is  easy  to  check  from  the  above  that  any 
symmetric  A  with  dim  ker  A  c  I  lies  on  the  SO (5) -orbit  of  a 

type  0,1,  or  2  matrix.  From  now  on  we  shall  call  these  orbits,  or 

any  representat i yes  of  them,  type  0,1,  or  2 . 

Finally  we  turn  to  the  remaining  A's  with  an  axis  of  equilibrium 
that  is  not  already  on  an  orbit  of  type  0,1,  or  2. 

We  say  A  is  type  5  if  dim  ker  A  =  2  ar.d 
f  A  =  0. 

•  I  f  A  type  5,  then  C  r-  sym  coir^jjjts  a 


6.10.  Definition. 


Proof.  S  =  QA  6  sym  implies  S'  =  A'  and  so  again  S  =  jf\, 
as  in  0.6  even  though  possibly  .V  V  =  0  .  In  this  case  Q  could 
be  any  rotation  about  il(..x).B 

All  the  foregoing  information  can  be  summarized  as  follows: 

6.12.  Theorem.  The  SO (3)  orbits  in  M_  fall  into  five  distinct 
types  according  to  the  way  in  which  they  meet  svm  jsec  Table  i  beiow.i  . 
Furthermore,  for  A  €  sym,  S^  =  • Q ' QA  '=  svm}  consists  of  1  {Q,} 
for  all  l  invariant  under  A  land  hence  S  '  sym)  except  in  the 
cases  (1)  dim  ker  A  =  2;  then  S^  also  contains  the  rotations 
through  any  angle  about  the  eigen- axis  of  A  corresponding  to  the 
non -zero  eigenvalue; 

(2)  A  =  0;  then  =  S0(3)  .  See  Table  2  below. 


Remarks 

1.  Table  1  highlights  the  fact  that  having  an  axis  of  equilibrium 

or  not  is  not  an  invariant  of  the  S0(3)  action  on  L  .  This  means 
there  are  equilibrated  loads  having  an  axis  of  equilibrium,  but 
which, when  rotated  globally  by  a  certain  amount  to  another 
equilibrated  load.no  longer  have  one. 

2.  Thus,  by  Theorem  5.1,  we  get  existence  of  solutions  to  the  traction 
problem  for  all  t ype s  of  astatic  loads  except  3,4.* 

3.  The  notion  of  tvpe  can  be  nulled  back  from  M_  to  L  with  a 

j 

little  care,  as  we  see  below. 


*In  particular  the  occurrence  of  0  solutions  by  Stonpelli  in  ~yoo  1 
is  seen  to  be  due  to  a  neglect  of  the  full  rotation  group  » see 
Section  S) .  Our  results  are  also  consistent  with  those  of  Ball 


35 


6 . 13.  Pet*  init ion . 

By  analogy  with  our  definition 

S  =  (Q  e  so (3)  !  QA  €  sym;  ,  A  €  M 
which  we  applied  when  A  is  of  type  0  ,  let  us  now  write 
S,  =  (Q£  SO (3)  ;  0^  €  Lg'  ,  i  €  L  . 

From  equivariance  of  k  we  clearly  have: 

6.14.  Lemma.  ^  . 

Note  that  the  map  -*■  0^  sym:  Ot~*  QA  is  an  embedding  for 
types  0,1,2  but  not  all  types  3,4,  because  of  the  isotropy. 
Pulling  back  to  L  ,  we  see  that  Qt—  01  is  an  embedding  S,  -  C ! 
if  k (Z)  is  of  type  0,1  or  2,  so  we  can  refer  to  2.  as  being  of 
type  0,1,  or  2  according  as  k(i)  is.  On  the  other  hand,  if  k(?.  1 
is  of  type  5  then 

either  (a)  0 ,  L  = 

or  (b)  C j  0  L  =  two  disjoint  circles  in  t  *  ker  k 

F:inallv,  if  k(-l)  is  of  type  4  then  C\  --  ker  k  3  L  and  any 

SO ( 3 )  orbit  C ,  is  all owab 1 e . 


ure  Set  eigenvalues 


ran  svc  r>a  1  intersection  <>t  t  lie  orbit  of  A  with  sym. 


types.  These  loads  are  all  pure  traction,  with  B  =  0. 


Tyne  1 .  Rotation  by  180° 
about  one  of  the  horizontal 
axes  produces  an  equilib¬ 
rated  load  with  no  axis  of 
equilibrium. 


Type  2 .  Any  horizontal  axis 
is  an  axis  of  equilibrium; 
vertical  axis  is  not  an  axis 
of  equilibrium.  Rotation  by 
180°  about  the  vertical  axis 
gives  an  equilibrated  load 
with  no  axis  of  equilibrium. 


Type  3  (cl)  .  The  load 
itself  admits  a  circle 
group  of  summetries 
about  the  axis  t-- 
which  is  thus  an  axis 
of  equilibrium. 
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Type  2  (  o) .  The  load  is 
not  symmetric,  but  the 
astatic  load  remains 
constant  under  rotation 
about  the  axis  — which 

is  thus  an  axis  of 
equilibrium. 


Tvpe  j  .  The  astatic 
load  is  zero:  all  axes 
are  axes  of  equilibrium. 


,oad  types 


'•  D 


3 T .  The  Bifurcation  Equation  and  its  Gradient  Character 

According  to  the  formulation  i  h  1]  of  our  problem,  we  wish  to 
solve  the  equation  H(\,A,n;Q)  =  0  for  Q,  where 

H f A ,A,n;Q)  =  SkewiQA)  -  \F ( \ ,Sym(QAl ,An ) , 


t.A.n)  is  near  (AA,n,,'l  <£  L  and  •  is  small.  Define  the  right - invari ant 
f  U‘  e  *■ 

vector  field  X  on  S0(3)  bv 
X0 

X  (Q.l  =  SkewiQA  }.Q; 

i.e.  right  translation  of  Skew(QAQ)  E  so(3)  =  TjS0(3)  to  T^S0(3) . 
Likewise,  we  shall  regard  H  as  a  right- invariant  vector  field  ci 
S0(3)  depending  on  the  parameters  ,\,A,n  by  setting 

Xp  ,A,n;Q)  =  HP  ,A,n  ;Q)  *Q 


Thus  , 


XfO.A(),n0;Q)  =  X  iQ). 

0 


Final lv,  note  that 


he  zero  set  of  \ 


3  ^  =  •  Q  S0(  31  Skew '‘./A,1  r- 

i' 


What  3 

'n 


for  various  t  vpcs  w.u  •; 


l  s 


■  i  .  on 


n  o  p  i 
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.  1 .  Lemma. 


Suppose 


sym  is  of  any  type.  Then  cor  0 


T  S 
Q  A . 


i.WQ|W  €  skew  and  WQA^  +  QA^W  =  0;  =  ker  DX^  (Q). 


Proof.  The  second  equalitv  is  clear  for  anv  A„ ,  because  r>XA  (Q)  : 

•  -  U  ;'n 

WQr—  skew  I'VQA  )*Q  .  For  the  first  one,  the  inclusion  C  immediately 

follows  by  differentiation  of  X  fQ)  =  0  in  0.  Equality  then 

A0 

follows  by  a  dimension  count;  recall  from  3.3  that  v  v  gives 

an  isomorphism  from  the  space  of  axes  of  equilibrium  for  A  (not 

necessarily  of  unit  length)  to  the  IV  €  skew  such  that  WA  +  A IV  =  0.  ■ 

T 

Recall  that  K  WQA^  +  (.QA  )  W  corresponds  to  the  linear  trans¬ 
formation  trace(OA^)  I-OA  ,  under  the  i  somoruh  i of  skew  =  SO  (31 

'O'  i) 

with  IR0.  When  Q  <=  S  ,  QA  is  symmetric,  so  this  t  rans  format  i  on 

A0 

is  symmetric  relative  to  the  Killing  form  on  SO ( 5 ) .  This  remark 
and  ~ . 1  proves  the  next  lemma. 


Lemma . 


.■suppose 


or  inv 


pe . 


Then  at  each  point 


of  S%  ,  the  ranvje  or  DX,  iQ):  T\S0(5i  —  T^SO'5)  is  the  orthoeonal 

_  Ao  - - -  \)  '5  Q  - - - 

complement  of  T . 


Next  we  recall  a  general  context  tor  the  bifurcation  of  vector 
fields  that  will  be  applied  to  our  situation  tcf.  Reeken  ’19"').  Let 
M  and  A  be  manifolds  and  X:  M  <  A  —  TM  a  smooth  vector  field  on 
M  depending  on  the  parameters  A  S  A.  We  seek  the  zeros  of  X.  For 
A  =  suppose  the  zero  set  S  of  X  is  a  known  smooth  compact 

submanifold  of  M.  Assume  that  M  carries  a  Riemannian  metric  and 
that  for  x  €  S,  the  range  of  D  X(x,.\q)  is  the  orthogonal  complement 
of  T^S.  The  normal  bundle  E  of  S  trivializes  a  neighborhood  U 
of  S.  For  each  x  €  U,  let  P  :  T  M  -*•  T  S  ,  ,  be  the  orthogonal 
projection  to  the  fiber  over  t(x)  ,  where  t:  E  -  S  is  the 

projection.  By  the  inverse  function  theorem,  there  is  a  unique 
section  :  S  -  E  such  that  P  X(-S,  (x)  ,.\)  =  0  for  x  €  S  and 
in  a  neighborhood  of  (assume,  for  example,  that  , '  is  compact  )• 

Let  Xfx,A)  be  the  orthogonal  projection  of  X(x.'-j  onto  the  tangent 
space  to  the  graph  of  i.  at  a  point  x  on  the  graph.  Thus,  XCx,'1 
is  a  vector  field  on  the  graph  of  f  and  finding  its  zeros  is  clearly 
equivalent  (for  small  '■l  to  finding  zeros  of  X.  We  call  the  equation 
X ( x ,  *  A  =  0  on  the  graph  of  •£.  the  bifurcation  equation.  Since  ^ 
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and  the  graph  of  are  diffeomorphic  under  ,  we  can  equally 

well  regard  X  as  a  vector  field  an  3.  This  reduction  of  the  nroolem 

is  often  known  as  the  Lianunov-Schnidt  method. 

The  above  procedure  may  be  applied  to  our  vector  field 

XI. \  ,A,n  :0.)  with  parameters  (\,A,n)  and  variable  x  =  ■}  €  SO (5)  =  M. 

Assume  \  is  near  cero  and  (A,n)  is  near  a  load  >'.\  ,nQ  •  where  A. 

is  of  arbitrary  type.  Thus,  there  is  a  unique  section  ,  of  the 

normal  bundle  to  3.  determined  by  the  Liapunov-Schmidt  procedure 
A0 

as  described  above.  Let  7(.\,A,nj  denote  the  graph  of  :. 

and  let  X(,\,A,n;Q)  be  the  orthogonal  projection  of  X  to  the 

tangent  space  of  7  at  Q.  Thus,  X  is  a  vector  field  on  7. 

As  above,  we  may  also  regard  X  as  a  vector  field  on  5.  . 

A0 

The  rest  of  this  section  is  devoted  to  proving  that  the  essential 
part  of  X  is  a  gradient.  In  the  general  context  above,  if  X  is 
a  gradient,  then  so  is  X  since  the  orthogonal  projection  of  a  gradient 
vector  field  to  a  submanifold  is  the  gradient  of  the  restriction.  This 
simple  version  does  not  directly  apply  to  our  situation  as  X  need 
not  be  a  gradient  vector  field  on  SO (35.  However,  the  "second  order" 
Taylor  approximation  X,  of  X  will  be. 

To  state  our  gradient  results,  recall  that  in  M  wo  defined  the 
quadratic  function  G:  —  skew  to  be  the  second  order  term  in  the 

Tavlor  expansion  of  F  about  0.  Thus  FIX,?)  =  V'i  •'  !  *  — G  (  ■  1  ... 

h 

where  G  is  a  quadratic  function  of  l.  The  appropriate  second  order 
approximation  to  the  vector  field  X  will  thus  be  defined  by 


— s  hew  ( 1 1  0 . 


x,r<  , A.ni'hl 


s k  ew ■ 0 A  1 
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Let  X,  be  the  second  order  approximation  of  the  vector  field  X  on 

S  obtained  by  the  Liapunov-Schmidt  orocedure.  Thus ,  is  the 

A0 


orthogonal  projection  of 


.3.  Theorem.  Suppose 


gradient  vector  field  on 


f(Q)  -  <  VQV 


X,,  onto  the  tangent  space  T  S  for  Q  *=  5,  . 

'  0  0 


Ajj  is  of  arbitrary  type.  Then  X.,  is  a 
S  .  In  fact,  X =  -grad  f,  where 


<v  -  <  v<V  *  t 


<7u0'c(7u0)>dV 


J8 

and  u„  =  Db(I  )  *(Q£-;i  i.e.,  is  the  uniaue  solution  in  C 

-  Q  <3'  '  0  -  Q  - - -  5>T7) 

of  the  linearized  equations  with  load  Qi,^  €  . 


Recall  that  the  pairing  between  loads  1  =  (B,t)  and  configurations 
(or  displacements)  is  given  by 

<:.,*)>  =  ^  B  (X)  •  C  (X)  dV  *  '  '  X  ■  :  X  dV  =  trace  k  ( 1 , 1 ) 

>  8  58 


and  physically  represents  a  potential  for  the  working  of  the  loads. 

T  t 

Observe  that  if  l  =  L  ,  then  1  ,0‘I.v  -  trace  ’  AO  i  -  tracefAO  )  =  \  '■  ,QI 

e  '  o 

for  all  Q  S  S0(3) . 

Remark .  In  the  second  term  of  X.,  and  f  we  can  replace 
bv  1,  .  However,  the  difference  is  higher  order,  so  •'  ,  is  sufficient 
for  subsequent  applications. 


Proof  of  7.3.  We  shall  show  that  X,  is  a  gradient  field  on 
S0(3)  which,  hv  the  remarks  follwoinc  ~.3.,  is  sufficient. 

We  proceed  in  two  parts.  Let  us  first  show  that  X  ‘Q1  is  the 
gradient  of  (  on  all  of  Sill 31  . 


TS 


7.4 

Lemma . 

Let 

£  L 

and 

A  =  k 

U).  . 

Define  a  vector  field 

on 

S0(3)  _bv 

X.iQ' 

—  Skew ( QA ) 

•Q  as 

above 

and  the  man 

of 

SO  (3)  _t£ 

A 

TR  bv 

■IUJ 

=  <  i 

,QTIo>  . 

.  Then 

X  =  -grad  I. 

Proof.  Two  simple,  but  useful  observations  are: 
If  B,  W  <E  M_ ,  with  W  S  skew,  then 


<  B , IV >  =  (  skew  B , IV  ) 


and 


If  B  €.  M  ,  l  <E  L  and  :  c  C  ,  then 


<  l,  B,i>  =  <  B,kU,j)  > 


To  prove  ".4,  we  compute  as  follows: 


dt(Q)*(WQ)  =  <  (WQ)TIg  ) 

=  <  O'.Oj1,  k(yg)>  by  (: 

=  <  (W0)T,.\  I 

=  < W  F , OA 


-  >V ,  skew  ; 0 \  i,'  bv  1  I 
=  WO,  skew  1 0-\)*  0  ■ 

=  - '  WO ,  \  ■ o 


T 
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This  deals  with  the  first  term  of  X.,.  To  deal  with  the  second 
term,  we  need  a  special  case  of  Betti's  reciprocity  theorem: 

7.5.  Lemma.  <  =  ^  <  WQ) , u^>  for  and  fWQ'l  2  £vm  . 

This  is  a  trivial  consequence  of  svmmetrv  of  D*(I„j  i.e.  of  the 

o 

elasticity  tensor.  It  is  also  proved  in  standard  references;  for 
example,  see  Truesdell  and  N’oll  ’  1965]  ;  p.  525. 


To  prove  '’.5.,  we  shall  also  need  to  calculate  the  second  derivative  of  the 

skew  component  of  i.e.  of  F(P)  =  Skew  f  k  ( J  (■;))]  .  Surprisingly, 

this  second  derivative  depends  only  on  the  classical  elasticity 

tensor  c  .  Recall  from  52  that  we  regard  c  as  a  linear  map  of 

1  T 

sym  to  itself  and  that  we  write  e  =  ^-(?u+(7u)  ). 

".6.  Lemma.  Let  F:  ”  -*■  skew  be  defined  by  Ff.f)  =  Skew  'k(vf'i))'. 

Then  F(Ig)  =  0,  DF(Ig)  =  0  and 


■>  f 

D“F(I  )(u,u)  =  2  Skew  (!  7u*cfe)dV  =  -2  Skew  k(?.  ,u) 
t-  JB  u 

where  =  (b^,T^),  b(l  =  -DIVfc(e))  and  7  =  c(e)*N  ,  Identifying 

skew  with  31'  ,  this  becomes 


-D“F(Ig) (u.u) 


b  x  u-  dV  +  f  t  -  u  dA 
u  J  u 

6  38 


IS 


Proof.  By  Lemma  4.2,  K  p)  =  Skew'1  PdV '  where  P  is  the 
/  ~ J  8 

first  Piola-Ki rchhof f  stress  tensor.  We  have  PfT  )  =  )t  so 

5 

f  .P 

F(I„)  =  0.  Also,  DF(I„)’u  =  skew  ,  4tr*7u  dY  =  skew  c*e  J\  - 
5  8  )gjF  ;R 

ti>  °  :  _ 

as  c*e  is  symmetric  and  since  —  (I_)  =  c.  To  compute  ii  - ,  »o 

jr  a 

shall  need  to  use  the  fact  that  S  is  symmetric.  Write  P  =  Ff  me 

use  the  product  rule  tc  obtain  DpP(F)*Tu  =  Tu * '3 i' F 1  +  FPCS( F;  •  .’u. 

Thus,  as  S(I  )  =  0, 
t> 


■> 

D'P  (  Ig)  •  (7u,7v)  =  7u*DfS^1s) -7v  *  7vDpS(FV7u  *  D'SUg)  * 1  7u,7v\ 


T  n 

Now  DrS(  Ip)  •  7u  =  D^S  (  L)  •  (7u  +  7u  )  =  c*e  and  D  S  {  I  n)  is  symmetric, 
r  b  ^  O  r  o 

so 


D"F(I„)«(u,u)  =  skew  !  d:T(I„)  (7u,~u)dV 
o  )  F  6 

8 


2  skew  (i  7u"ci.e)  JV 


8 


Finally,  this  equals 


-2  skew 


U.  dV 


&  \j>,  dA  } 
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by  the  divergence  theorem,  so  the  last  statement  follows.  T 


Example.  For  i  homogeneous  isotronic  material, 


:  i e  !  =  '  i trace  e ) I  -  2ue 


where  e  =  [7u  *  I7u)  j/2  and  a  are  the  Lame  moduli. 


D“F(I_) tu.u)  =  2  Skew  (  "\7u* f trace  ( .u) ] I  *  2u7u*e}dV 

6  3 


=  2  Skew  \  [  trace  (7u)  ]  7u  +  u7u*7u;dV  .* 

‘  8 


Let  us  next  see  what  -.b  says  about  the  quadratic  term  G  in 

the  Taylor  expansion  of  F.  For  t>  £  C  we  have  the  i dent i tv 
r  svm 


FU)  =  FP  *($)) 


where  Pg:  L  -*•  Lg  is  the  orthogonal  projection.  Thus,  as  OF  and 

OF  are  zero  at  I„  and  0  resnect  i  vel  <•,  and  P  LV'(ITj  =  Df(I_!  , 
b  e  5  S 


D"Fd_)fu,v)  =  D“F(0)(M(IJ*u,W(IJ.vl 

O  DO 


Hence  for  ;  £  and  writing  u,  =  D<M  T  1  we  have  the  identity 

e  .  8 ' 


G(i)  =  D“F(I_) (u,  ,u,1 
o  •.  < 
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-GU)  =  2  Skew  [  [  b  ®  u  dV  +  I  t  ®  u  dA] 

L  J  *  J 
S  38 


=  2  Skew  k((b,x),u  ) 


where  b  =  -DIV (c*  (el  , ) ,  and  T  =  c(e,)-N  and  e^  =  [7u,  (7u,)  ]  '2  . 

However,  these  last  equations  say  exactly  that  (b,i)  =  l ,  and  so  we  get 


-  jGU)  =  Skew  kU,u  ) 


Completion  of  the  proof  of  7.3.  The  derivative  of 
\  T 

Q  b-  (  Zq,  7Q  u^>  in  the  direction  WQ  is  given  by 


<  J.Q..  t(WQ)Tuq)  *  (  ;Q,  ^Tuwo  > 


<  ( ((VO)  uQ )  (by  Betti  reciprocity,  ".51 


•  <  W  ,k(Qf,0,uQ)> 


by  f 2) 


■  <  V!  ,  Skew  k(Q?.p,u^)>  by  (1) 


-<  WO  ,  Skew  k(Q;o,u(.j)Q  ) 


(  WO,  >*0  > 


by  i  5 )  | 


2 8 .  Bifurcation  Analysis  for  Type  I 

We  now  discuss  the  solutions  of  the  basic  equation 

H( a , A ,N ;01  =  Skew  (OA)  -  \F  (a  ,  Sym'OAl  ,Qn)  =  0  1  1 

for  the  load  i  =  (A,n)  near  a  load  =  (A  having  an  axis  of 

equilibrium  and  of  type  1,  and  for  A  near  n  .  We  shall  also  obtain 

the  stability  of  the  solutions  and  finally  we  shall  compare  our  results 

with  those  of  Stonpelli  [1958].  For  type  1  we  need  to  do  a  bifurcation 

analysis  on  the  circle  S  corresponding  to  the  degenerate  zero  set 

A0 

of  H  when  A  =  0  and  1=1^.  The  analysis  has  some  features  in 

common  with  the  papers  of  Hale  [1977]  and  Hale  and  Taboas  [1981], 

Without  loss  of  generality  we  can  assume  that  A^  =  diap ( i ,  -  a , -c ) 

where  0  t  a”  t  .  Thus,  from  §6  the  set  5.  of  zeros  of  Skew  ('OA')  for 

A0 

Q  £  SO f 3)  is  given  explicitly  by  the  following  two  points  and  circle: 


S  =  r di agf 1 . - 1 , - 1)  diagf-1 ,1 ,-l) }  U  C 

An  0 


Cl 


where 


r  '■  11  \ 

q  -  I  y  x  u  J  I x  =  cos  A ,  y  =  s in 
'()  0  1  / 


The  loads  corresponding  to  the  two  points  are  A  =  Jiag  i  a  ,  a ,  c  and 
+■ 

A  =  diag i -a , -a  ,c  )  . 

From  ".5,  we  are  led  to  study  the  critical  points  of 
T  T 

f(Q)  =  (  ,  ,i)  ,  n,0  >  on  .  Note  from  the  divergence 


theorem 


hat 


Thus  the 


<  v°  uo }  =  ;p'7uo*  c'eo}  >dv 

where  u  =  D$(I)  1  (02.^)  and  =  [.u^  *  fTu^,)‘]/2  .  Thus  the 

function  f  is  computable  from  linearited  elasticity  alone,  which 

leads  to  the  curious  observation  that  our  "second  order"  nonlinear 

elasticity  here  involves  no  more  data  than  linear  elasticity,  but 

merely  processes  the  information  in  a  different  way.  Writing  Q= 
x  -y  0  \ 

y  x  0  ]  as  in  (2),  f  becomes  a  nolvnomial  of  degree  2  in 

v  0  0  1  J 


1  oy 


(x,y).  Write  the  two  terms  of  f  as 


f(Q)  =  =  lh0  +  V  +  b;y)  +  |Carx-  *  a2xy  *  a.y-  *  ajX  *  a.v  - 


f  (x.y)  =  ff(.x,y) 


[4) 


which  defines  the  numerical  constants  and  . ab- 

define  new  parameters  by  writing 


and  letting 


(S' 


f  f x ,y)  =  3jX"  *  a,xy  ♦  u..y-  *  u,x  *  a-y  * 


Note  that  x  . . .  .  ,ntfi  depend  on  our  parameters 
the  elastic  moduli  of  the  material.  Thus , 


as  well  as  on 


L\  =  3r  a:  =  az ’  ‘s  =  a3 


=  ^1  *  V  l3  =  ^2  *  V  S'  5  rb''  ‘  •'* 


Replacing  ^  by  Q£  ,  where  Q  is  as  in  2',,  effects  a 

rotation  of  the  x-y  plane.  Thus,  by  rotation  of  if  necessary, 

we  can  assume  x.,  =  0  . 

Let  us  fix  :t.  ,  t.  and  consider  the  bifurcations  of  ceros  of 


dr* 

-r—  -  Z  ( X,  -X.  XV  * 

dr  0  1  • 

i-x 

ci  ,. 

-  a ,  v-  on  L>  i.e., 

,  critical  points  of 

o 

4'  ’ 

1 

on  5*1  with  a. 

and 

a_  as  parameters. 

o 

Set  M  a  f  fa  ,a_  ,■']  €  R"  <  S*  ^4—  =  O'-  the  manifold 

4  o  dr  4  ^ 

of  critical  points  of  f * .  Indeed,  M  is  a  manifol  i  and  can  be 

parametrized  bv  ,z:  R  <  S1  M  ,  pfu.f)  =  (-2l'a  +u)  cos  -,  -2(i.*uiS  ^  ■ 

lorn 

">  i  2 

Denote  by  m:  1R  x  S  -*■  IR  the  projection  onto  the  first  tactor. 

■>  ->  i 

3.1.  Lemma.  Set  l  =  f2(a1-a.)“  -  a  -  108a  “a_“  ( a  -a.'i"  .  f~  j 

-  -  1  o  as  4slj> 

-> 

If  a,  -  'i.  i  0  ,  then  t  .  M  -<■  1R"  is  a  proner  stable  nan  in  I  i , ,  a.  1 -snace , 

-  I  J  -  - * - * - 4 -  4  j  - 

and  its  set  of  critical  values  is  the  astroid  defined  by  .  =  0 
i  see  Figure  below  )  . 


df* 
d  ■ 


Since  the  number  of  points  in  ”  ft)  (i.e.,  the  terns  of 

at  a  =  (  t,,a_)'l  is  a  constant  over  '  ■  n  or  •  0  ,  we 
4  j 


obtain 


S.2.  forollarv. 


df* 


has  4  zeros  if 


•  ;) 


and  has  2  zero; 


0 


Proof  ot'  Lemma  3.1. 
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The  critical  set 
3(4,9)  SR  <  31  |  ot1 

critical  values  of  m 


Z  of  t->o:  R  <  S2  -*•  R“ 

7  o 

sin~9  +  a,  cos~3  +  y  =  0}  . 
can  be  parametrized  by 


is 

Thus,  the  set  of 


a ,  =  -Z ( a , -a,)  cos0  9 
4  1  o 

a_  =  2(ctj-a_)  sin0  9 


(8) 


Since  I  consists  of  4  cusp  points,  4  fold  lines  and 
T’C;!  is  1-1  ,  by  a  result  of  I'Tiitney  (see  Mather  [1969]  or 
'^olubitsky  and  :uillemin  [1973]),  one  knows  that  Ttop  is  a  stable  map. 

Eliminating  9  produces  the  bifurcation  set 

C(a  -a  )]2/3  =  a  2/3  -  u,2/3  (9) 

1  ^  4  D 


For  -  ou  i  0  ,  (8)  describes  the  astroid  shown  in  Figure 


F igure  4- 


Next ,  observe  that  for  real  numbers  A,  fi  and  C, 


A  +  B  +  <J  =  0  if  and  only  if  A'1  +  R'1  *  (V5  =  3ABC 


by  virtue  of  the  identity  A'5  *  R'1  ♦  C"'  -  a  ARC  =  ■  A  ♦  R  »  Ci'A” 

AB  -  BC  -  CA)  .  Applying  !  101  to  f'A)  shows  that  ,  i equivalon 

:  ;  :  2/3  2/3 

~i .  *  t.  -  -la,  -  x.  )  =  -art,  x.  -  i_!l 

A  .1  l  a  A  .>  1  a 

gives  the  stated  conclusion.  B 


Cub i nr  both 


The  family 


of  functions  on  S  with  parameters 


*  j  > 

4  o 


so 


-if* 

d6 


eniovs  a  universal  property.  Consider  a  perturbed  family 

df*  -  m 

-tt—  +  g(.\,c,e)  ,  with  g(0,0,-5''  =  0  tor  (i,cj  -  IR  <  R  .  To 

do 

,  jf* 

each  (\,c),  denote  by  M,  _  =  \  (j5~  *  g )  !  \ ,  c ,  ,  x.  ,t) :  = 

the  "manifold"  of  zeros. 


S.3.  Lemma .  For  (\,c)  sufficiently  small,  the  sets  M.  _  are 

manifolds  and  there  exist  two  smooth  families  of  diffeomorohisns 

a,  :  R-  ~  R~  and  4\  :  M,  -  M  defined  for  \,c  sufficientlv 

■  ,  c  -  \  ,c  '  ,  c  -  - — 

small,  such  that  nf.  =  •>,  n  ,  and  T  =  identitv,  vy  =  identity. 
- -  a  ,  c  •  X ,  c  -  0,0  ■  0 , 0 


1  ’ 

Proof:  For  \,c  sufficiently  small,  the  map  t IR  •  S  *  M.  . 

V  ,  c  '  ,  C 

o,  (y.-a)  =  (-2(a1+y)  cos  ft  +  sin  f  4|-  ,  -2f:t.  +  u)  sin  J  -  cos  •-  ,  -  ■  I 

defines  a  parametrizat  ion  of  .  3y  Lemma  S.l,  'r’C.  R  <  -*  R 

is  an  unfolding  of  the  proner  stable  man  .  Thus,  one  can  find 

_ !  2 

diffeomorphisms  ?*.  ^  on  R  «  S‘  and  R”  respect  i  vely 

such  that  •;«-  f-r'c,  )  =  ’f*.  •  This  lemma  follows  bv  letting 

s,c  '  ,  c  '  ,  c 

i>,  =  ,:c  ?*,  ..-.j1  .  | 

•. ,  C  *  ,  C  '  ,  c 


how,  we  are  ready  to  state  our  main  result  on  the  number  if 

solutions  near  type  !  loads.  Let  ■  =  ■' 1  c  ■  ieperui  smoot a  1  y  -mi  i 

parameter  in  Rm  ,  with  "  -  =  ■  Recall  that  d  ; - 

iet'ined  bv  ~  •  .  a,  ,  a,  and  a.  by  1  i  . 

*  > 


3.4.  Theorem .  Let  oe  a  type  1  load  with  k  i  ;  =  .  a,-a,-c!, 

-»  •> 

0  i  a"  i  c”  ,  a.  =  0  ,  and  a,  i  a_  .  Then,  there  exists  a 
_  -  1  o  - 

(smooth!  function  1(a,c1,  A(a,0)  =  A(a^,a_!  +  0('-J  defined  tor 

I \ , c 1  sufficiently  small  and  \  >  0  such  that,  the  traction 
problem  has: 

( i  1  four  solutions  for  the  load  \ .i ( c )  i_f  A ("•. , c !  0 

i two  of  them  near  C  )  . 

An 

(ii)  six  solutions  for  the  load  \A(c)  if.  A(.l  ,c)  >  0 

(two  of  them  near  Cl  . 

A0 

Proof:  The  bifurcation  of  zeros  of  a  fcf.  3~)  on  C.  (=S1j 

A0 

is  the  same  as  finding  zeros  of  ■ -X,  4-r  >  =  —  •+  —  g(\,c,f) 

=  or  ^  +  g  ,  where  g(0,0,i)  =  0  .  Let  r 

atr.  -  “  at*  * ,  c 

be  the  family  of  diffeomorphisms  found  in  Lemma  S.3.  Take 
b  (4(c)!  b,(>.(c)) 

A(''. cl  =  AC1.  ,  (  — - -  +  a  ,  -  *  a.!  which 

*•  9  C  '  4  —  '  J 

has  the  desired  property.  1 

Next,  we  want  to  determine  the  "generic"  structure  of  the  bifurcation 

set  K  =  -A  =  0:  in  (A,c)  space,  '•  >  0. 

For  m  =  0  ,  A(ot,,a„)  =  A  (a,,  a.)  *  0  .  it  is  clear  that  K  -  0  . 

4  o  4  a 

Indeed,  our  traction  problem  has  two  solutions  near  C  it 

Ao 

AC  a  .a.  1  0  ,  and  four  solutions  near  C  if  !a  ,a„)  l»  . 

1  j  \3  *  -1 

b  Afc',  b.Vcl 

F:or  m  =  1  ,  consider  k:  c  i — *! — ^ -  *  a  ,  — — r -  +  a.) 

where  \|c  1  is  the  linear  part  of  ;.!c  ’  .  Th  i  <  represents  i  straight 
line  which  we  assume  to  intersect  the  astro  id  transversely  if  they 


) 


meet.  \'otice  that  K.  -  -o'  <  T*. .  c  >*=  •<: 


.5  tne  inverse  image  ot 


the  astroid  : defined  bv  eanation  (9)),  under  the  man 
b^lc)  '  '  b,l(cl 

h,  :  c • — ii.  (— =r -  +  a  ,  — a.  l  .  I’j,  as  in  the 

A  ‘  y  C  -A  4  D  A  ,  C 


proof  of  Theorem  3.4).  Write  \c  =  c  .  Recall  that  4(c)  =  A(cj  +  0(  c.~) 


and  consider  the  map 


b .  A  ( c ) 

c  K-n,  (,\,c)  =  ^ - 

'  .  V  C  - 


0(X), 


b-,A(c) 


Of-.))  . 


Since  the  astroid  is  bounded  and  "j.  is  close  to  the  identity, 

\,c 


there  exists  an  interval  such  that  K, 


'Ac  S  K.  }  C.  i-M,M) 


for  \  >  0  ,  c  sufficiently  small.  Applying  the  isotopy  theorem 
for  transversal  maps  (see  e.g.,  Hirsch  (1976))  to  the  family  h- 
through  Iiq  =  k  ,  we  conclude  that  the  bifurcation  set  ~X  consists 
of  0,  2,  or  4  curves  with  slopes  k  ^(astroid)  (see  Figures  £,&)  . 
Thus,  for  example,  by  choosing  c  i  0  sufficiently  small,  and 
letting  \  -*■  0  (i.e.,  consider  the  load  a4(c)),  one  can  pass  from 

a  parameter  region  where  there  are  two  solutions  near  the  circle 
(four  in  all)  to  one  where  there  are  four  near  the  circle  (six  in 
all)  and  back  again  to  the  two-solution  region  isee  Figure  i ). 

Such  a  situation  is  not  dealt  with  in  the  analysis  of  Stoppelli  [1958], 


For  m  >  2 


let  us  suppose  that  the  affine  nap: 


b.A(c)  bnA(c) 

c*-  1 — - -  +•  a  ,  -  +  a.)  is  onto,  where  again  Afc) 

—  4  „  D 

is  linear  part  of  tie)  .  Without  loss  of  generality,  we  may  also 
assume  that  b^A(c)  =  c ^  and  b_,A(c)  =  c.,  ,  where  c  =(Cj,c,,a'l  . 


Notice  that  K ,  _  =  {  (Cj,c,)  (\,Cj  ,c,,;)  €  K  ;  is  the  inverse  image 


of  the  astroid  under  the  map  h^  .. 


h1  M  ■’  ) 


bp  CO 


+  a.  =  a. 


„  <  —  +  a  ,  +  ar  )  .  Set 

'  *  ^  t  *  L  -i  *  *  “  *+  -  V  d> 


•  .  .\c„  -  c_  .  and  consider  the  rr.ap 


As  be tore,  \> 


\,=  =  '  tci,c:)  '  C1'-'C: 


is  bounded  unitorml 


for  \  >  0  ,  c  sufficiently  small.  Xpplying  the  isotopv  theorem 
for  transversal  maps  to  the  family  h,  _  through  the  affine 

isomornhism  h,,  .  we  conclude  that  the  bifurcation  set  is  a 

0, 0 

cy 1 inder- 1  ike  set  along  the  i-axis  with  base  a  cone  over  the 


astroid  in  Cj.c,  space.  The  first  order  approximation  of  this 
cone  is  given  by  the  cone  over  the  astroid  in  the  plane  \  =  1, 
centered  at  f-2a,,  -2a.  1  with  ’’sire’'  4. a, -a.  (see  Figure  7). 

4  j  1  j 


7 


Figure  7 


Next  we  discuss  the  stability  of  the  solutions  correspondin';  to 


loads  near  load  of  type  1.  This  can  be  determined  by  combining 
our  stability  results  for  type  )  Theorem  5.5.;  together  with 
well-known  stability  results  for  the  cusp.  We  make  the  same 
assumnt ions  as  those  in  Theorem  3.4. 


3.5.  Theorem.  Let  =  diag(a, -a, -c) ,  A^  =  diag(a,a,c.l  and 
★  ■* 

Aq  =  diag(-a, -a,c)  as  above .  The  indices  of  the  bifurcating 

solutions  are  boxed  in  Table  5.  iRecall  that  stable  solutions  have 

index  =0.)  In  each  case  the  circle  represents  C  defined  by 

0 


equation  (2) . 


Note  that  two  stable  solutions  bifurcate 
c  >  a'  .  In  all  other  cases  the  solutions 
unstable . 


off  the  circle  when 
near  the  circle  are 


L 
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which  is  a  constant  (independent  of  x,  y ).  In  this  situation, 

=  :t  =  0  and  so  our  theorems  do  not  apply. 

3.5.  Example .  Consider  the  same  traction  problem  as  above,  but 
with  a  homogeneous  nonisotropic  hyperelastic  material  whose 
linearized  elasticity  tensor  is  given  by  c(<-)  =  £  -  ^  Jiag  £  . 


In  this  case, 

u  (X)  =  c'1  Fqf 0 

0 

-a 

0\-i 

o'  x  , 

,.fy 

-v  O' 

x  o  ; 

Q  L_  \  0 

0 

-c  /- 

l  0 

0  1  / 

where  c  ^ (F)  =  F  +  diag  F  .  Then 

<  V  QV  =  j^7vc(7V>  dv 

/ Tax  ay  0  \  /ax  ay  0  \  \ 

ay  -2ax  0  )  ,  [  av  -ax  0  i  y 

Vo’  0  -2c  ;  \  0  0  -c  / 

7  ->  7  ^  -) 

=  4a"x"  +  2a“y"  +  2c“ 

1  ■>  .  ,  , 
Hence,  \  -  3a  >  0  ,  and  our  traction  problem  tn  ■  ■  nas 

six  solutions  ('four  near  )  ,  with  stability  determined  by¬ 


table  5. 


Next  we  shall  discuss  how  to  obtain  the  results  of  stoppel 1 i 
1958  .  as  a  special  case  of  our  analysis,  We  refer  the  reader  to 
the  statements  of  Stoppelli's  results  in  Grioli  19t>2,  p.  58..  In 
this  approach  one  focuses  attention  on  bifurcations  that  occur  on  the 
circle  by  examining  what  happens  at  a  particular  location  on  the 
circle  and  •  «  ■. n  .  V.’e  can  assume  this  'joint  is  il,9)  i.e.  r  -  ) 
with  no  loss  of  generality. 

First  of  all,  if  .t-,  *  u-  ?  0  then  (1,0)  is  not  a  critical 
point  of  f * ,  so  there  are  no  solutions  near  (1,0).  We  may  assume 
then  that  un  +  u.  =  0  and  then  the  Taylor  expansion  of  f*  about 
9=0  becomes 


* 

o  1 


(higher  order 


of  Grioli  19<i.2  . 


■] 

i 

Case  2.  If  - 

+  ;t- 

-  -7-  =  '>  and  X-, 

f  u,  then 

df  * 

w  = 

5  .  2 

+  0(-?\i  and  so  : 

here  are  ■  r  -  -u 

ut  i  >)ns 

i fold  point ] . 

This  is  Theorem 

I  on  p.  3>  of  >r:oi 1 

i.u'2 

T 

Case  5.  If  -•< 

3 

-  '4  ■ 1;  -  0 

1 1  . 

but  t.  -  *  -r—  *  1 

1  .1  0 

! 

.  d  f  *  4 

taen  ITT  =  IU1 

*”3 

-  iii-' . 

so  there  are  4  or  . 

solutions  !  cusp  point.).  This  is  Theorem  H  on  p .  SS  of  irioli  T.b-J 


Furthermore,  it  we  express  our  constants  x.  t-a.  1  :n  terms  of  t in¬ 
elasticity  tensor  c  and  solutions  of  the  linearized  problem 
using  (3)  above,  we  find  the  same  conditions  for  these  three  case  ■ 
as  is  given  on  p.  5"  of  Grioli  I'.lhJ.. 

Thus  we  recover  the  results  of  Stoppelli  on  type  1  loads.  Vs 
was  explained  in  the  introduction,  however,  this  analysis  is  only 
local  on  the  circle  and  does  not  give  the  full  story  of  the 
bifurcation  picture,  even  in  this  case.  The  compl<  .e  hi  fmvat  i  a 
analysis,  including  stabi  1  i  ty .  is  summarized  by  our  Figure  .  uv. 

Table  3. 
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